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Two-dimensional cluster-cluster aggregation is studied when clusters move both diffusively and 
sediment with a size dependent velocity. Sedimentation breaks the rotational symmetry and the 
ensuing clusters are not self-similar fractals: the mean cluster width perpendicular to the field 
direction grows faster than the height. The mean width exhibits power-law scaling with respect to 
the cluster size, {r x ) ~ sS L = 0.61 ± 0.01, but the mean height does not. The clusters tend to 
become elongated in the sedimentation direction and the ratio of the single particle sedimentation 
velocity to single particle diffusivity controls the degree of orientation. These results are obtained 
using a simulation method, which becomes the more efficient the larger the moving clusters are. 
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I. INTRODUCTION 

Aggregation of particles and particle clusters is still 
of great interest not only due to the number of ap- 
plications it has in chemical engineering, material sci- 
ences, and atmosphere research but also due to its fun- 
damental role as a simple model system for growth un- 
der non-equilibrium conditions 0, 111- The effects 
of the interplay of diffusive and ballistic motion in the 
case of fractal aggregates has been considered only re- 
cently [lllElalIli|lS|ll|ll|ll. The relative 
strength of these two mechanisms will vary with cluster 
size. This is the case, for example, in colloidal suspen- 
sions, where both the diffusivity and the sedimentation 
velocity are usually expected to depend algebraically on 
cluster size. In this paper we study the effect of these 
two processes on the cluster structure. 

The present understanding of colloidal aggregation 
under gravitation is as follows. According to experi- 
ments jg, sedimenting clusters do not rotate, their 
anisotropy is independent of their size, and they have 
no preferred orientation. Based on a scaling relationship 
between cluster size and sedimentation velocity the sed- 
imenting clusters are argued to be self-similar and their 
fractal dimension to be significantly larger than that in 
diffusion-limited aggregation. The restructuring of parti- 
cles inside a cluster caused by hydrodynamic stresses was 
claimed to be the reason for sedimenti ng c lusters being 
more compact. Recent simulations 0, [H^ have shown 
that even without restructuring it is possible to produce 
clusters with an apparent fractal dimension that is close 
to the value observed in experiments, if there is a ve- 
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locity difference between clusters of different sizes. The 
cross-over to the sedimentation-dominated regime can be 
measured by directly observing the average size of the 
aggregates; however in this regime the fractal dimension 
still follows DLCA-like scaling p| . 

In this article we study scaling properties of clus- 
ters, which are formed in cluster-cluster aggregation 
when both the diffusivity and the sedimentation veloc- 
ity depend algebraically on cluster size. This induces a 
cross-over from diffusion-limited aggregation to a process 
where large clusters grow while settling by aggregating 
smaller clusters from a time-dependent size distribution, 
n s (t). We are mainly interested in generic features and 
hence concentrate on two dimensions, in which the effects 
of sedimentation on cluster structure can be studied with 
the most ease numerically. To compare our results to the 
existing ones, we also consider other structural character- 
istics on a qualitative level. For quantitative studies one 
should study more elaborate models with hydrodynam- 
ics and other delicate issues involved. The main focus 
here is on the self-similarity of aggregates, which is the 
crucial assumption made in the analysis of data in pre- 
vious, similar studies. For this purpose, we consider the 
scaling of four radii as a function of cluster size. These 
are indicated in figure ^ 

Our main result is that sedimenting clusters are not 
self-similar fractals. Hence, an algebraic relationship - 
as assumed a priori here - between cluster velocity and 
its mass does not necessarily imply self-similarity of the 
aggregates. The mean cluster width, considered in the 
direction perpendicular to the sedimentation velocity, 
grows algebraically as a function of cluster size but the 
mean cluster height does not follow a simple power-law, 
at least not in the size range considered. It is an open 
question as to whether the asymptotic scaling of the clus- 
ters transverse to the velocity follows the mean cluster 
width, or, vice versa; much larger simulations are needed 
to settle this issue. Also the principal radii of gyration, 
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FIG. 1: Illustration of the principal radii of gyration r m i n and 
TWx calculated from the radius of gyration matrix T and the 
radii r x and r y , which point perpendicular and parallel to the 
sedimentation velocity (shown by the thick arrow), respec- 
tively. The longer side of the rectangle shown equals to the 
maximal distance between any of the particles of the cluster 
and the angles f> and 6 are used to measure the orientation 
with respect to the direction of the sedimentation velocity. 



given by the eigenvalues of the radius of gyration matrix 
(see Fig. ^ and Eq. (J5J), grow with rates that are not 
equal. We present a scaling argument that explains how 
the width grows algebraically but simultaneously leaves 
the origin of the actual value of the exponent open. The 
fact that this particular quantity exhibits scaling brings 
up some interesting issues, like where does the universal, 
Peclet-number independent value originate from? 

The difference in the growth rates indicates that clus- 
ters can not be described simply using the fractal di- 
mension and methods relying on isotropic scaling prop- 
erties. For example, studies using the radius of gyration 
give only an effective measure for the cluster structure. 
Anisotropic growth also implies that the cluster shape 
and its structure change in time, and we consider the 
former in particular in detail. We also analyze the orien- 
tation of clusters and discuss the similarities and differ- 
ences between our results and previous experimental and 
numerical studies of the same phenomenon 0, 0, IToj . 

This paper is organized as follows. The model is de- 
fined in section [H] A simulation method becoming the 
efficient the larger clusters become is introduced in sec- 
tion [ffi] The scaling properties using the radius of gyra- 
tion matrix are studied in section ITVl Section [V] consid- 
ers the anisotropy and orientation of clusters. Section lVTI 



concludes the paper. 



II. MODEL 

Consider fractal clusters aggregating in a suspension. 
The sedimentation caused by gravitation may be ig- 
nored for aggregates with characteristic radii smaller 
than about 1 /im. In this region the essential features 
of the growth are well described by the diffusion-limited 
cluster-cluster aggregation (DLCA) model The 
idea is that clusters move diffusively, and the diffusion 
constant depends algebraically on cluster size: D(s) ~ 
s 7 . For a cluster diffusing in a quiescent fluid it can be 
argued that the diffusion exponent 7 = — 1/df, w here df 
is the fractal dimension of the cluster [iiJ 113 ■ When- 
ever clusters (particles are clusters of size one) collide, 
they irreversibly aggregate together. No cluster restruc- 
turing is allowed. 

Both the dynamics and the structure of clusters formed 
in the DLCA are well understood 0, ^| . Defining the 
number of clusters of size s at time t as n s (t), the 
cluster size distribution obeys dynamic scaling n s (t) — 
S(t)~ 2 f(s/S(t)), where the exponent —2 follows from 
mass conservation and the average cluster size S(t) = 
J2 S s2n s/^2 s sris ~ t Z w ith z being the dynamic expo- 
nent. In two dimensions the fractal dimension of clusters 
is d f = 1.44 ± 0.02 0[l9| and including rotational dif- 
fusion or varying the value of 7 has no essential effect on 
it P, ^| . The anisotropy measured using the ratio of the 
principal axis of gyration is about 2.4 [20j. 

Here the DLCA model is considered in the presence 
of the sedimentation of clusters. The force exerted on 
a cluster sedimenting in a fluid consists of the buoyancy 
force Fb = VApg and the viscous drag force Fd = —Cv, 
where V is the volume of the object (for an aggregate 
consisting of s particles it is s times the volume of a 
single particle), Ap = p p — pf is the difference between 
particle and fluid densities, g is the acceleration due to 
gravity, v is the sedimentation velocity and C is a posi- 
tive constant depending on the specific form of the clus- 
ter [23. For a ball of radius r the constant C = 6irr]r, 
where rj is the kinematic viscosity. For complicated, say 
fractal-like objects, C is generally unknown and we make 
the usual assumption that the aggregate will behave as 
a compact object with an effective hydrodynamic radius 
Rh, i.e. C ~ Rh im I23. If one further assumes 
a scaling relation between Rh and s, the sedimentation 
velocity depends algebraically on cluster size v{s) ~ s s , 
where 5 is the sedimentation exponent. For a fractal one 
would have s ~ Rfj resulting in S = 1 — 1/df but here 
we just take this as an assumption since it is not a priori 
clear that the ensuing clusters are fractals with respect 
to drag resistance. As 7 < and S > aggregation 
is diffusion-limited for small clusters but becomes domi- 
nated by sedimentation for large ones. 

The model described above neglects many issues re- 
lated to real suspensions such as hydrodynamic interac- 
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tions, restructuring and rotation of clusters by diffusion 
or at aggregation to mention a few. However, our pur- 
pose is not to try to simulate all aspects of sedimenta- 
tion driven aggregation but rather elucidate what kind 
of universal behavior one could expect to have. Hence, 
the results reported should be considered on a qualita- 
tive level when comparing to experiments. We discuss 
further these issues and the assumptions behind 7 and S 
in the conclusions, in the light of the simulation results 
obtained. 

A useful measure for the relative strengths of diffusion 
and sedimentation is the Peclet number |23| 

P c (s) = v(s)r/D(s), (1) 

where r is the radius of a single particle and v = \v\. As 
the Peclet number depends on cluster size, in the follow- 
ing we use the one particle Peclet number Pe = f c (l). 
Figure|21shows examples of clusters formed in simulations 
for different values of Pc. 



III. COMPUTATIONAL ASPECTS 

The simulations are done on a two-dimensional lat- 
tice with periodic boundary conditions. Initially a con- 
centration 4> of lattice sites are filled randomly. One 
filled lattice site is considered as a single particle. In 
the dynamics a particle or a cluster consisting of par- 
ticles is selected randomly and time is incremented by 
Nft)^ 1 ^^, where N(t) is the number of clusters at 
time t and £l max is the maximum mobility of any of the 
clusters at that time. The cluster mobility is defined as 
^( s ) =Pu( s )+Pit( s ) + 2 P-L( s ); where Pu(s), pj.f(s), and 
p±(s) are proportional to the probabilities for a cluster 
of size s to move in the direction of the field, opposite to 
it and to the directions perpendicular to it, respectively. 
The selected cluster (of size s) is moved if a; < fi(s)/ ^max, 
where x £ (0, 1) is a random number selected from a uni- 
form distribution. If the selected cluster is not moved the 
simulation proceeds by selecting another cluster, other- 
wise another random number y £ (0, 1) is selected from 
a uniform distribution. The cluster is moved to the field 
direction if y < pu(s)/n(s), opposite to the field direc- 
tion if p u (s)/n(s) < y < \pu(s) +p iT (s)]/Q(s), and 
otherwise to one of the directions perpendicular to the 
field with equal probabilities (— p±(s)/n(s)). If after the 
move the cluster has no overlap with others the move is 
accepted and another cluster is selected randomly. Oth- 
erwise the move is taken back and the cluster is aggre- 
gated with all the clusters it attempted to overlap before 
selecting a new cluster for the next attempt. 

The terms Pn(s), p|f(s), and p±(s) are of the form 

p u (s) = (4C^ + C lS 6 + Cf S 25 )/2 

p u (s) = (AC^-C lS s + Cy s )/2 (2) 

p±{a) = 2C ± s\ 



where C|, Cm, and C± are non-negative constants chosen 
such that Pii(s) > for all cluster sizes. This choice gives 
for the sedimentation velocity v s — C^s 5 and for the dif- 
fusivities Du(s) — C^s 1 and D±(s) = Cj_s 7 , in the direc- 
tions parallel and perpendicular to the field, respectively. 
In the following we consider only the case Cn = C± and 
denote D(s) = D^(s) = D±(s). Physically, only the ratio 
of the sedimentation velocity to the diffusion constant is 
relevant and hence we report the results as a function of 
the one particle Peclet number Pc = Cj,/Cj_. The ex- 
ponents are taken to be 7 = —0.70 and 5 — 0.30, which 
are obtained from relations (5=1+7 an d 7 = — 
using the fractal dimension of DLCA clusters df = 1.44. 
This is of course a valid choice initially, when DLCA is 
the dominating process. It should be underlined that the 
main point in these values is that eventually sedimenta- 
tion will prevail, and that the both transfer mechanisms 
obey an algebraic dependence on the cluster size s. 

It is known that in DLCA cluster structure, and es- 
pecially the effectuve fractal dimension, depends on con- 
centration 28] . In simulations reported here, the concen- 
tration is kept fixed to (f> = 0.01. The DLCA calculations 
are done on a lattice of size L11 x L± = 4000 x 4000 
and when sedimentation is included 15000 x 2000 unless 
stated otherwise. To minimize finite size effects, all the 
simulations are run under the conditions that none of the 
particles has traveled the size of the lattice in the sedi- 
mentation direction and the size of the largest clusters is 
smaller than half of the lattice size in any direction. 

The above choice for Pu(s), Pii(s), and p±(s) assumes 
that sedimentation does not have any effect on the dif- 
fusion properties of a cluster. This may not be the case 
as the velocity fluctuations of sedimenting, compact, and 
non-aggregating particles are larger in the sedimentation 
direction than perpendicular to it pil I2M [2(J . The de- 
pendence of the strength of fluctuations on aggregate size 
is unknown, so we make the choice of isotropic fluctua- 
tions. The precise form of fluctuations should be unim- 
portant as the upper critical dimension for sedimentation 
driven aggregation is expected to be one [27J . 

Next we shortly discuss the implementation of the al- 
gorithm. It is constructed to be efficient for low concen- 
trations and for large cluster sizes. Although the method 
is described using a two-dimensional lattice, the general- 
ization to higher dimensions and the continuum case is 
straightforward. 

Each cluster is represented by a virtual rectangle, 
whose size is the smallest possible one to enclose the 
cluster such that two sides of the rectangle are paral- 
lel to the sedimentation velocity (see Fig. |2J). The sites 
belonging to a cluster are stored with respect to the, say, 
lower left corner of the enclosing rectangle, for two rea- 
sons. The first one is to minimize the time needed to 
move a cluster. Regardless of the cluster size a move can 
be performed by updating one single value i.e. either the 
x- or y-coordinate of the lower left corner of the enclosing 
rectangle. The second reason deals with the collisions of 
clusters. It is fast first to check if the rectangles overlap 



FIG. 2: Typical clusters formed in simulations for Pe = 0.01, 1, and 10 (from left to right). They have 41167, 87585, and 80512 
particles and are contained in rectangles of size 688 x 942, 1169 x 1307 and 1373 x 1625, respectively. The sedimentation 
velocity points downwards. 



each other and only if they do to check for the overlap of 
the clusters. 

Obviously it would be very time-consuming to com- 
pare the rectangle that is moved to all the other rect- 
angles. Hence, the whole system is divided into rect- 
angular blocks and clusters are identified to belong into 
blocks according to the rectangles enclosing them. Then 
only clusters belonging to the same block(s) have to be 
checked. 

Figure |31 compares CPU-times using the implementa- 
tion described above to that with a "traditional" algo- 
rithm, where clusters are represented as occupied sites 
on a lattice and moved site by site. For average cluster 
size smaller than about ten the traditional algorithm is 
faster as one needs only to compare the occupancy of 
a few sites. On the other hand, for large clusters the 
method using enclosing rectangles becomes much faster, 
being independent of cluster size, and usually only a cou- 
ple of comparisons are needed to check the overlap of a 
cluster with others. In the "traditional" implementation 
one needs of order S(t) comparisons, which slows down 
the computation considerably when S(t) becomes large. 

In the simulations reported here the "traditional" 
method is used for S(t) < 10 and the "rectangle" one 
for S(t) > 10. Hence, most of the computation time 
is used when the average cluster size is small. To ob- 
tain good statistics from the interesting, large cluster 
size region, we further apply a cloning method. It has 
proved to be efficient in many applications including, for 
example, studies of percolation clusters, native states of 
polymers, and reaction-diffusion systems |29j |. The ba- 
sic idea is simple: at a fixed average cluster size S cop y 
one makes n cop y copies of the cluster configuration and 



continues the Monte Carlo simulation independently for 
each of these systems. We use S copy = n copy — 10. We 
tested that the results obtained using copying were in- 
distinguishable from the ones obtained without it. The 
averages reported in this article are taken over 40 runs 
without copying and 10 runs which were copied 10 times. 
To increase the scaling window further would necessitate 
enormous CPU resources; even with the algorithm pre- 
sented above we estimate that to augment the maximum 
cluster size by one order of magnitude would require of 
the order of one CPU-year on a fast computer. 

Note that in the range of s- values shown, the typical 
scaling of S(t) is not algebraic manifesting a cross-over 
from diffusive to sedimentation-dominated dynamics |27| . 
Since in the asymptotic regime the mean-field dynamic 
exponent z = 1/(1 — 5 — l/d e g), for any reasonable value 
of d e ff (the effective fractal dimension) z will have a value 
much larger than unity, and the cross-over effects are 
noticeable since the DLCA value is much smaller. 



IV. SCALING PROPERTIES OF CLUSTERS 

Complex aggregates are often characterized be their 
fractal dimension df, which can be found by the relation 
of the cluster size to its characteristic radius 

s~{R) d ', (3) 

where the brackets denote averaging over clusters of size 
s. Usually one considers the radius of gyration 
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FIG. 3: CPU-time as a function of the average cluster size 
for 'traditional' (O) and 'rectangle'-algorithms (A). Each 
symbol represents the time used from the previous inspection 
time. The inset shows the ratio of the two CPU-times. 



where denotes the position of the ith particle and 
f CM = i Ylt=i i s the center of mass of the cluster. 

As the field breaks the rotational symmetry it is not 
obvious that clusters will scale isotropically. Hence, we 
consider the radius of gyration matrix 

\ ± yx ± yy J 

where 

1 s 

T a/ 3 = - y^(n iCC - rcM,a)(n,p - r C M.p) (6) 
s f-f 

%— 1 

and x and y refer to the x- and y-components of the 
position vectors. The larger and smaller eigenvalues of 
T, denoted by A max and A m i n , respectively, are related 
to the radius of gyration by Rg = A max + A m i n - The 
square roots of the eigenvalues, called the principal radii 
of gyration, r maX:inin = A^f x min , may be considered as 
lengths of an ellipsoid describing the cluster shape and 
the corresponding normalized eigenvectors e max and e m i n 
define the directions of the principal axes. The preferred 
direction is given by the eigenvector corresponding to the 
larger eigenvalue. We further consider the radii in the 
directions perpendicular and parallel to the field, which 

1/2 

are related to the T by the relations r x = T xx and r y — 
Ty!y 2 \ respectively, and to the radius of gyration by Rg = 
r x +r y . Figure^depicts the four radii for a cluster formed 
with Pe = 0.1 and of size 90110. 

Figures|U(a) and 0(a) show the growth of the average 
radii defined above as a function of the cluster size for 
Pc = 1 and Pe = 10, respectively. Although the growth 
seems rather similar for different radii, the local slopes 
(shown in Figures 4 (b) and 5 (b)) change continuously 
with the cluster size. Only the radius perpendicular to 
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FIG. 4: (a) The mean radii (r max > (O), (r mi n) (□), (r*) (♦), 
and (r y ) (T) as a function of the cluster size for Pe = 1. (b) 
The local slopes. 



the field obeys nice scaling of form (r x ) ~ s' x with l x — 
0.61 ±0.01. 

The overall behavior for other Peclet values (Pe = 0.01 
and 0.1; not shown) is rather similar and, especially, the 
value of l x is independent of the Peclet number. For 
comparison, Figure shows the scaling of the radius of 
gyration [Eq. J3J] as function of cluster size for various 
Peclet numbers. The local slopes reveal that only for 
the DLCA there is a scaling relation between these two. 
For a non-zero Peclet number one has to consider the 
scaling using the radius of gyration matrix [Eq. or 
some other observable not relaying on isotropic scaling. 
One should emphasize that for all Pe > any asymptotic 
behavior is still far away in the range of accessible radii 
of gyration; the actual values reflect more the cross-over 
in the cluster shape and internal structure (density) than 
a tendency to e.g. become compact in the limit t — » oo. 
This can be underlined by comparing with the scaling of 
the quantities, in Figs. 0] and [SJ 

The reasons for the existence or lack of scaling with 
respect to r x and r y can be considered by a simple scaling 
argument. If r x scales with s, then also its derivative with 
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FIG. 5: (a) The mean radii (r max > (O), <r min > (□), (r x ) (♦), 
and (r y ) (T) as a function of the cluster size for Pe = 10. (b) 
The local slopes. 




FIG. 6: Cluster size as a function of the radius of gyration 
for the DLCA (O), Pe = 10 (■), Pe = 1 (0) ; Pe = 0.1 (*), 
and Pe = 0.01 (A). The inset shows the corresponding local 
slopes indicating the change in the effective fractal dimension 
of sedimenting clusters. 



respect to s does so. We have that dr x /ds « Ar x /As 
and 

Ar x =P(r x )x Sr x (7) 

where P(r x ) denotes the probability that r x changes and 
Sr x measures the typical change in r x per such an aggre- 
gation event. 

Consider a cluster of mass si, which collides with an- 
other one, with S2 < s\. Then Sr x ~ r x {s2) ~ Sj? 
and P(r x ) ~ r x (s 2 ) /r x (si) ~ (s 2 /si)' x . Thus dr x /ds — 
s% /(si S2) an< i by making the ansatz {.S2) ~ asi (since 
S2 is cut-off by si, a < 1), we obtain that 

drjds ~ s'f- 1 , (8) 

exactly as one should if the growth is algebraic. 

In the case of r v , the probability P(r y ) does not have 
to be algebraic. This is simply so because (see Fig.0) the 
roughness of the lower part of a cluster growing via sedi- 
mentation will also play a role: many of the smaller clus- 
ters will be "swallowed" inside the fjords opening down- 
wards. It seems feasible that this leads to a change in 
the compactness of aggregates, which in turn results in 
a non-algebraic relationship between r y and s. This is 
in fact a "restructuring" process, related to the changing 
compactness of aggregates. 

Note the analogy with single cluster growth in various 
ballistic growth models (with shadowing, see Ch. 5.7.2 
of 30]). There it is known that the bulk of the cluster 
becomes compact (df = d), and that the boundary of 
the cluster undergoes interesting roughening behavior. In 
our case, l x does not show signs of such a cross-over (to 
compact geometry) though l y might do so asymptotically. 

The equivalence to cluster growth is by no means clear 
in our case since the "deposited" clusters are sampled 
from the n s (t) at each t, the distribution of which is 
broad (Fig.0demonstrates an example, see also ref. |3l|). 
These clusters have a typical width-to-height ratio de- 
pending on s, making the "deposition problem" much 
more complicated than such ones usually are. This will 
remain true for any t, as the n s (t) remains non-trivial. 
The implication is that for any finite Peclet number there 
should be a slow cross-over process that depends self- 
consistently on the n s (t) and on the cluster structure 
that has been formed through sampling smaller clusters 
from the same. In any case, the fact follows that the 
cluster structure is anisotropic: the more recently aggre- 
gated parts of clusters are not scale-invariant with the 
parts that have been created earlier in the aggregation 
process. 

Consider next the local exponents. For typical cluster 
sizes the local exponents of (r ro i n (s)) and (r x (s)) seem to 
be close to the same value (see Figs. 0] (b) and|5] (b)). 
Also the slopes of (r max (s)} and (r y (s)) behave similarly 
for large cluster sizes. This indicates that the cluster are 
elongated in the field direction, i.e., (r y ) > (r x ). The 
elongation is more pronounced for high Peclet numbers. 
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FIG. 7: The averaged (40 runs) cluster size distribution n a (t) 
for Pe = 1, at t = 1024. 



The same conclusion is obtained by directly consider- 
ing the angular distribution using the eigenvector corre- 
sponding to the larger eigenvalue (see Section IVjl . 

Although in the size region studied (r y ) > (r x ), the lo- 
cal slope of (r x ) is larger than that of (r y ) indicating that 
clusters grow faster in the direction perpendicular to the 
field. This implies that very large clusters should become 
elongated perpendicular to the sedimentation direction. 
Rough estimates for the crossover values are s = 8000, 
80000, and 300000 for Pe = 0.1, 1, and 10, respectively, 
and the corresponding radii of gyration 80, 200 and 500. 
These values are, unfortunately, larger than we are able 
to simulate but should be achievable in experiments, see 
for example 8J. The most important conclusion, is how- 
ever a choice between two asymptotic scenarios. The 
data implies that asymptotically (r x (s)) will scale with 
the maximum radius. Hence, should this be the case, 
there are two possiblities to begin with: either the max- 
imum radius also starts to follow the (r x (s)}, with its 
scaling exponent (0.61 above), or then the asymptotic 
scaling of both could become volume-like (~ s ' 5 ). 



V. ANISOTROPY AND ORIENTATION 

Next we consider how the cluster shape evolves with 
time. To have an easily comparable measure for clus- 
ters of different sizes, we consider the scaled half- width 
of a cluster as a function of its scaled height Y = 
(y - ymin)/(y m ax - 2/min), where y max and y min are the 
maximum and minimum y-coordinates of any of the par- 
ticles in the cluster. The half-width at height Y is defined 
as 

W(Y) = l ^xOO-ZminC^ (Q) 
2 Umax 2/min 

where x max and x min are the maximum and minimum 
^-coordinates of any of the particles at that height. 
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FIG. 8: Average cluster half- width as a function of scaled 
height for Pe = 0.01 (open symbols) and 10 (filled symbols). 
The clusters are divided into four size categories: 2500 — 
5000 (O), 5000 - 7500 (□), 7500-10000(0), 10000 - . . . (A). 



The relative width of a cluster increases with size as 
can be seen from Figure [HI in which the half-width is 
shown for two Peclet numbers. Note, that the relative 
width may well become larger than the height of a clus- 
ter. For Pe = 0.01 the clusters are symmetric with re- 
spect to the point Y = 0.5 but for Pc = 10 they are wider 
at the bottom. This is due to the fact that as large clus- 
ters sediment faster than small ones they gather mass at 
bottom side. In other words, the upper part of a cluster 
is shielded by its lower part. For large clusters this results 
in a triangular-like shape, a particularly nice example of 
which is shown in Fig. [21 The sedimentation driven ag- 
gregates are sparse and have complicated scaling prop- 
erties (see Figs. 0] and QjJ) as the clusters formed in the 
diffusion-dominated regime are themselves fractals. It is 
an open question as to what is the best characteristics of 
this kind of anisotropic cluster shapes. One possibility 
would be to look at the lower and upper parts separately 
(e.g. the widths thereof, after splitting the cluster at the 
center of mass). 

The orientation of sedimenting clusters is considered 
using the angle tp = arccos(e y ■ e max ) between the eigen- 
vector corresponding to the larger eigenvalue and the 
unit vector in the direction of the external field. Fig- 
utcEJ shows the angular distribution for two values of the 
Peclet number. 

To study the effect of cluster size on orientation the 
clusters are divided in five different size classes. There is 
no difference in orientation with respect to cluster size. 
The main point, however, is that the clusters are ori- 
entated such that they prefer to have the longer princi- 
pal axis aligned with the field and that the orientation 
distribution depends on the Peclet number. For exam- 
ple, for Pe = 10, there are practically no clusters with 
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FIG. 9: Angular distribution with respect to the sedimenta- 
tion direction for Pe = 1 (open symbols) and Pe = 10 (filled 
symbols) and for several size classes: 1000 — 2500 (O)j 2500- 
5000 (□), 5000-7500 (0), 7500 - 10000 (A) and 10000 ... 
(V). 



ip > 45°. This is not in contrast to the conclusion made 
using the half-width as there is an important difference 
between these two approaches, when one averages them 
over cluster sizes. The averaging of the principal radii 
gives information about the typical cluster shape whereas 
averaging cluster width takes into account the selection 
of a preferred direction. We believe, that if one could 
simulate larger cluster sizes, another peak would start to 
grow at <p = 90°. 

To consider our simulations in the light of the 
three-dimensional experiments we also calculate cluster 
anisotropy and orientation as in Ref. 0. First, a clus- 
ter is enclosed in a rectangle, whose longer sides give the 
maximal distance between any of the points in the cluster 
and the shorter ones just touch the cluster (see Fig. QJ. 
The cluster "radius" is defined as R = (L + Z)/4, where L 
(I) denotes the length of the longer (shorter) edge. The 
orientation is measured using the angle 9 between the 
sedimentation velocity and the longer edge. 

Figure ^| shows the average anisotropy and the aver- 
age orientation angle as a function of R. The averages 
are taken over intervals of length one, which explains the 
noisiness of the data for large cluster sizes. There is a 
clear transient region, where both quantities vary with 
cluster size but for large clusters they become constant 
within statistical errors and the value of the anisotropy 
ratio L/l ~ 1.4. This is the same as the value of three- 
dimensional clusters [8j. However, since r x and r y have 
different scaling properties this would be a bit surpris- 
ing, if it were the asymptotic value. It seems most likely 
that the constant value reflects more the fact that the 
statistics becomes sparse, with only a few clusters in this 
range of R. One can compare the anisotropy to the ratio 
of r max and r m i n ; this starts for s small from the DCLA 
value and is thus different. Another explanation might be 




300 400 
R=(L+l)/4 



500 
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FIG. 10: a) Average orientation angle (8) and b) aspect ratio 
(L/l) as a function of radius R = (L + Z)/4 for Pe = 1. Solid 
lines represent sliding averages using 30 data points. 



the crossover from field-elongated aggregates to ones that 
are wider in the direction perpendicular to the sedimen- 
tation velocity: At R — 300 the cluster size s is approx- 
imately 45000, which is of the same order of magnitude 
than the value estimated for the crossover (r x ) w (r y ) in 
Sec. ITVl (.s=80000). 



VI. CONCLUSIONS 

We have studied the structure of two-dimensional clus- 
ters, which aggregate at contact and whose motion is 
dominated by diffusion for small cluster sizes and by sed- 
imentation for large ones. Considering the eigenvalues 
of the radius of gyration matrix, we show that due to 
the sedimentation the aggregates grow slower in the di- 
rection of the sedimentation velocity than perpendicular 
to it. Moreover, only the cluster width scales with clus- 
ter size and hence the clusters are not self-affine either, 
at the very least within the size range currently accessi- 
ble in simulations. We underline that the data implies 
a cross-over scale, at which the maximal cluster radius 
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will begin to scale with the average cluster width, but no 
further conclusions can be drawn about this regime. 

Interestingly, the scaling exponent characterizing the 
growth of the width is independent of the Peclet num- 
ber, which characterizes the strength of the influence of 
sedimentation on aggregation. It is the one that defines 
"self-similarity" in the growing aggregates, but only in 
a limited sense since it concerns only a one-dimensional 
characteristics of the clusters (in contrast to e.g. a frac- 
tal dimension). The width scaling argument presented 
leaves open what the value of the width exponent would 
be but hints about the reasons as to why scaling is ob- 
tained in this respect, only. Here the clear conclusion 
is that since the Peclet number has no influence on the 
value (s=s 0.61) it appears to be a universal one, and as 
such should be measured experimentally. The value to 
be obtained from three-dimensional experiments would 
hence be independent of the viscosity of the suspending 
liquid but not necessarily take the same value as obtained 
here. Recent experiments have shown that the radius of 
gyration can be used to distinguish between the DLCA 
and sedimentation regimes, with a clear change in the 
rate of increase with time However, the fractal di- 
mension was shown to initially remain at the correspond- 
ing DLCA value. We believe this underlines the fact that 
in such simulations as ours it is easier to concentrate on 
measures that characterize the behavior of the large ag- 
gregates, as for instance R g (s). 

Let us next discuss the justification of the simulation 
rules with the asymmetric growth in mind. As far as 
the growth is dominated by diffusion, the aggregates 
will be self-similar and diffusion characterized through 
D(s) ~ s 7 is reasonable. When sedimentation starts 
to dominate this no longer holds but then the precise 
characterization of diffusion is unimportant anyway. In 
retrospect the algebraic dependence of the sedimentation 
velocity on the cluster size is just a lucky choice, a poste- 
riori justified by the fact that r x scales algebraically. On 
the other hand, the relation S = 1 — 1/ dt is invalid. This 
raises the academic question that how does the value of 
the universal exponent l x depend on the values of 7 and 
81 One may argue as follows. As the fractal dimension 



of DLCA clusters is independent of 7, it probably does 
not change the scaling properties of sedimenting clusters 
either. However, 6 may and probably will have an effect; 
at least on the universality grounds there is no reason 
why the value of l x would not change with S. It would of 
course be possible to attempt much more elaborate and 
time-consuming simulations, in which the hydrodynamic 
radius of each aggregate is computed self-consistently. 

The anisotropic growth is also seen when studying the 
average form of aggregates. Larger clusters are always 
relatively wider than their smaller counterparts. For 
small Peclet numbers the clusters are rather elliptic but 
they became wider at the "bottom part" due to a shield- 
ing effect for large Peclet numbers. We also studied the 
orientation of clusters with respect to the direction of 
the sedimentation velocity. Clusters are oriented such 
that their height is larger than the width, but eventu- 
ally the situation should be the opposite as the cluster 
width grows faster than the height. The larger the Peclet 
number the stronger the cluster orientation anisotropy. 

The inequality of the growth rates in different direc- 
tion implies that the scaling properties of sedimenting 
clusters can not be studied using the radius of gyration, 
which implicitly assumes the clusters to scale similarly 
in different directions. As there exist data from experi- 
ments and simulations 0, Il0j | considering aggregation 
of sedimenting clusters, it would be worthwhile to check 
for the possibility of anisotropic growth in these cases, 
too. Experiments or simulations with larger aggregates 
could also reveal if the height of sedimenting aggregates 
will eventually scale with their size. It would definitely 
be also of interest to make simulations in the three di- 
mensional case, and if cluster restructuring or rotational 
diffusion is allowed. 
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